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Heat  transfer  of  phase  change  material  (PCM)  in  an  open  cell  micro-foam  structure  was  numerically 
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micro-foam  is  a  body-centered-cubic  (BCC)  lattice  embedded  with  spherical  micro-pores.  Two  different 
simulation  methodologies  were  applied.  One  is  the  high-fidelity  direct  numerical  simulation  (DNS), 
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PCM/micro-foam  systems,  with  the  effective  thermal  conductivity  derived  from  direct  simulations  and 
expressed  as  a  power  law  of  porosity. 
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1.  Introduction 

Metal  or  graphite  foams  [1]  filled  with  phase  change  materials 
(PCM)  are  attractive  for  thermal  energy  storage  (TES)  [2]  particu¬ 
larly  in  high  heat  flux  applications  as  they  provide  both  high  stor¬ 
age  capacity  and  short  time  constants.  The  latent  heat  of  PCM 
provides  high  heat  capacity  whereas  the  foam  structure  increases 
the  effective  thermal  conductivity.  In  earlier  studies  of  phase 
change  materials  in  foams,  the  volume-averaged  (macroscopic)  ap¬ 
proach  such  as  one-temperature  and  two-temperature  models  has 
been  used,  in  which  either  local  thermal  equilibrium  (LTE)  be¬ 
tween  PCM  and  foams  was  assumed,  or  the  simplified  Newton’s 
cooling  law  was  applied  [3,4].  The  volume-averaged  simulation 
can  be  effective  for  porous  structure  with  simple  geometry,  how¬ 
ever,  this  method  can  be  restrictive  for  the  modeling  of  foams  with 
complex  geometry.  The  complex  foam  structure  which  affects  the 
effective  thermal  conductivity  and  interstitial  heat  transfer  coeffi¬ 
cient  and  subsequently  the  PCM  melting/solidification  needs  to  be 
accurately  modeled  to  provide  a  better  understanding  of  the  heat 
transfer  behavior  of  the  system.  The  conventional  method  of  effec¬ 
tive  thermal  conductivity  derived  from  packed  bed  with  discrete 
solid  particles  as  is  done  for  simple  porous  structures  cannot  be 
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applied  here  without  modification  and  more  accurate  methods 
such  as  direct  numerical  simulation  are  required. 

In  this  paper,  we  investigate  the  phase  change  in  a  low  porosity 
micro-foam  using  both  volume-averaged  and  direct  numerical 
simulation  [5].  The  foam  was  modeled  as  an  ideal  structure  with 
BCC  unit  lattice  and  spherical  micro-pores.  The  temperature  at 
the  top  surface  of  the  porous  structure  was  kept  constant.  The 
objective  of  this  work  is  to  compare  various  simulation  methodol¬ 
ogies  and  identify  one  that  can  model  the  heat  transfer  process  in 
PCM  and  micro-foam  structure  with  sufficient  efficiency  and  accu¬ 
racy.  The  volume-averaged  simulations  were  carried  out  using 
one-  and  two-temperature  models,  whereas  in  the  DNS,  separate 
energy  equations  were  applied  to  micro-foam  and  PCM  domains 
and  conjugate  heat  transfer  was  considered  on  the  interface.  Re¬ 
sults  including  melted  PCM  volume  fraction,  temperature  differ¬ 
ence  between  PCM  and  micro-foam,  instantaneous  temperature 
profile,  and  heat  flux  are  discussed. 

2.  Model  and  methods 

2.1.  Physical  model 

The  physical  model  and  geometry  of  the  simulated  micro-foam 
is  shown  in  Fig.  1.  A  cube  with  the  side  length  of  2.56  cm  was  used 
for  volume-averaged  simulations  (Fig.  la).  The  top  surface  was 
held  at  several  constant  temperatures  (Ttop  =  334,  351,  or  368  K), 
and  the  other  surfaces  were  insulated.  In  direct  numerical  simula¬ 
tion  (Fig.  lb),  the  micro-foam  geometry  was  composed  of  a  series 
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Nomenclature 

Alphabetic  symbols 

q" 

heat  flux 

Ajs 

area  of  PCM/micro-foam  interface 

Greek  symbols 

Cp 

specific  heat 

a 

power  factor  in  Progelhoff  model 

H 

side  length  of  cube 

P 

density 

hfs 

interstitial  heat  transfer  coefficient 

y 

porosity  (0-1) 

hy 

volumetric  heat  transfer  coefficient 

ST 

temperature  difference  between  micro-foam  and  PCM 

k 

thermal  conductivity 

Subscripts 

L 

latent  heat  of  fusion 

eff 

effective  value 

rd 

distance  from  the  cell  centroid  to  the  top  surface 

Ini 

initial  condition 

s 

ratio  of  PCM  surface  area  to  the  total  area  on  the  top 

Top 

top  surface  of  the  simulation  domain 

boundary 

s 

solid  (micro-foam) 

T 

temperature 

f 

fluid  (melted  PCM) 

Tm 

melting  temperature  of  PCM 

t 

time 

of  unit  cells,  and  each  was  represented  by  a  solid  cube  (side  length 
of  400  jim)  with  nine  spherical  pores  arranged  in  a  BCC  lattice  [5- 
7].  The  length  in  Y-direction  was  kept  the  same  as  that  in  the  vol¬ 
ume-averaged  simulations  and  comprises  of  64  unit  cells.  The  X  &  Z 
lengths  are  much  smaller  (0.02  cm)  to  keep  the  mesh  size  and  com¬ 
putation  time  at  an  affordable  level.  The  top  surface  of  the  unit-cell 
array  is  specified  with  the  same  temperature  condition  as  that  in 
the  volume-averaged  simulations.  The  bottom  surface  of  the 
unit-cell  array  was  insulated.  Symmetric  conditions  were  specified 
at  other  surfaces.  The  initial  temperature  of  PCM/micro-foam  sys¬ 
tem  was  300  K  in  both  volume-averaged  and  direct  simulations. 

We  considered  aluminum  foam  and  docosane  (paraffin  wax 
with  C22H46  composition)  PCM  whose  melting  temperature  of 
317  K  is  below  the  top-wall  temperature.  It  was  assumed  that 
the  micro-foam  is  fully  saturated  with  docosane  and  there  is  no 
volume  change  of  PCM.  All  thermo-physical  parameters  of  these 
two  materials  (Table  1)  were  kept  constant  in  the  simulations. 

The  foam  porosity  was  0.757  and  pore  size  is  360  pm  for  the 
current  simulations.  The  fluid  flow  within  the  melted  PCM  region 
can  be  neglected  because  such  low  porosity  micro-foam  provide 
a  dominant  resistance  force  to  suppress  the  flow  driven  by  the  nat¬ 
ural  convection  [8].  Thus,  only  the  energy  equations  were  solved  in 
the  simulations. 

2.2.  Direct  numerical  simulation  of  PCM  in  foam 


Table  1 

Thermo-physical  parameters  for  PCM  and  micro-foam. 


P  (kg/m3) 

cp  (kJ/kg-K) 

k  (W/m-K) 

L  (kj/kg) 

Tm  (K) 

PCM 

785 

2.89 

0.4 

260 

317 

Micro-foam 

2719 

8.71 

202.4 

N/A 

N/A 

2? 

Qd  S3 

II 

<1 

<1 

1 

d) 

dTs  w 

PsCP,s  dt  ~  V  ' 

<1 

(2) 

Conjugate  heat  transfer  was  considered  at  the  PCM/micro-foam 
interface: 


Tf  =  Ts 
kf§='< 


(3) 


Subscripts  /  and  s  in  Eqs.  (l)-(3)  stand  for  PCM  and  foam  respec¬ 
tively;  p,  cp,  T,  t,  k,  and  L  are  density,  specific  heat,  temperature, 
time,  thermal  conductivity  and  latent  heat  of  fusion.  In  the  right 
hand  side  of  Eq.  (1),  f  is  the  volume  fraction  with  the  value  of  1 
in  the  melted  PCM  region,  0  in  the  solid  PCM  region,  and  0-1  in 
the  mushy  region.  Effects  of  porosity  and  pore  size  have  been  di¬ 
rectly  included  in  the  ideal  foam  geometry. 


Without  the  need  for  extra  ad  hoc  assumptions  such  as  local 
thermal  equilibrium  (LTE),  direct  numerical  simulations  are  con¬ 
sidered  as  much  closer  to  the  real  situation.  In  DNS,  various  energy 
equations  were  applied  to  PCM  and  the  foam  domains  respectively, 
as  shown  in  the  following: 


2.3.  One -temperature  simulation  of  phase  change  in  foam 

One-temperature  model  assumes  that  PCM  and  the  micro-foam 
have  the  same  temperature.  Therefore  only  one  energy  equation  is 
applied  for  the  whole  system: 


Fig.  1.  (a)  Geometry  and  thermal  boundary  conditions  in  volume-averaged  simulation;  (b)  ideal  foam  microstructure  (PCM  not  shown)  in  direct  simulation.  The  top  surface  is 
held  at  a  constant  temperature. 

2 
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(yPfCpj  +  (1  -  y)PsCp,s )  -ft  =  v  •  (keffVT)  -  pfyL (4) 

where  y  is  the  porosity,  with  value  1  standing  for  pure  solid.  The 
effective  thermal  conductivity  keff  is  a  function  of  porosity,  thermal 
conductivity  of  each  component,  and  more  importantly,  the  detail 
structure  of  the  foam.  The  simplest  and  commonly  used  expression 
kejf  =  ykf  +  (1  -  y)ks  is  essentially  the  arithmetic  mean.  In 
PCM/micro-foam  system,  a  more  accurate  keff  can  be  obtained  by 
experiments  [9-11],  or  by  DNS  of  real  or  ideal  foam  geometries 
[5-7,12-14].  In  this  paper,  we  have  carried  out  DNS  to  find  the 
effective  thermal  conductivity,  which  matches  the  Progelhoff  model 
[15]  keff  =  ykf  +  (1  -  y)aks  with  a  =  1.31  through  curve  fitting,  when 
pores  in  ideal  foam  structure  is  arranged  in  a  BCC  lattice  [16]. 

2.4.  Two-temperature  simulation  of  phase  change  in  foam 

In  two-temperature  model,  energy  equations  for  PCM  and  foam 
need  to  be  considered  separately  because  of  the  local  thermal  non¬ 
equilibrium  (LTNE).  Their  heat  transfer  is  modeled  by  the  Newton’s 
law  of  cooling: 

ypfrs 5  =  V  ■  ( W VTf)  -  Pf)’L |  -  hv(Tf  -  Is)  (5) 

r)T 

(i  -  y)PsCP,  =  v  •  (fcs4rvrs)  -  hv(T,  -  Tf)  (6) 

where  effective  thermal  conductivities  for  PCM  kf>eff  and  for  the 
foam  kSteff  are  still  related  to  foam  structure.  Here  we  used  two  sets 
of  effective  thermal  conductivities:  one  corresponds  to  the 
arithmetic  mean  such  that  kf:6ff  =  ykf  and  ks^f  =  (1  -  y)ks.  The  other 
corresponds  to  the  Progelhoff  model,  with  kffig  =  ykf  and 
ks,eff  =  (1  -  y)aks  respectively.  The  volumetric  heat  transfer  coeffi¬ 
cient  hv  =  AfshfS ,  where  the  specific  surface  area  Afs  is  the  ratio  of 
interstitial  surface  area  to  the  representative  elementary  volume 
(REV),  and  h/s  is  the  interstitial  heat  transfer  coefficient  between 
PCM  and  foam.  Both  hfs  and  Afs  are  highly  dependent  on  the  foam 
structure.  In  this  work,  the  specific  surface  area  A/s  =  1.08  x  104  m_1 
is  obtained  from  the  foam  structure  described  in  Section  2.1.  The 
interstitial  heat  transfer  coefficient  is  defined  as  hfs  =  q'fS/(Ts  -  Tf), 
where  q'f  is  the  heat  flux  at  the  fluid/foam  interface,  Tsis  the 
volume-weighted  average  temperature  for  micro-foam  and  Tf is 
the  mass-weighted  average  temperature  for  fluid.  The  interstitial 
heat  transfer  coefficient  hfs  is  calculated  in  a  similar  manner  as  that 
by  Krishnan  et  al.  [6,7].  We  applied  direct  simulation  on  the  unit¬ 
cell  level,  and  based  on  the  calculated  heat  flux  across  the  interface 
and  the  average  PCM  and  foam  temperatures  we  obtained  the  value 
of  interstitial  heat  transfer  coefficient  hfS  =  9.09  x  103  W/(m2  K). 

2.5.  Problem  setup 

For  these  simulations,  the  commercial  CFD  package  ANSYS 
FLUENT  13.0  was  used  with  double  precision.  Parallel  and  batch 
processing  were  carried  in  the  AFRL  DoD  Supercomputer  Resource 
Center  (DSRC).  To  apply  the  two-temperature  model  in  FLUENT,  a 
user-defined  scalar  (UDS)  was  used  for  micro-foam  temperature. 
The  convergence  criterion  of  absolute  scaled  residual  for  energy 
equation  was  set  to  be  10-9  for  volume-averaged  simulations 
and  10-11  for  DNS.  Grid-  and  time-step-independent  results  were 
ensured  first.  We  found  that  the  grid  with  90  k  cells  and  time  step 
size  dt  =  0.02  s  are  enough  for  one-  and  two-temperature  simula¬ 
tions,  while  we  used  the  grid  with  9  million  cells  and  dt  =  0.008  s 
for  direct  simulations.  A  typical  DNS  needed  over  30,000  CPU  hours 
(i.e.,  192  CPU  cores  for  around  160  h),  much  more  expensive  than 
volume-averaged  simulations. 


3.  Results  and  discussion 

The  direct  numerical  simulations  show  very  complicated  micro¬ 
scale  evolutions  of  melting  front  and  temperature  distributions, 
which  are  very  different  from  the  one-dimensional  macro-scale 
behavior  distributions  generated  by  the  volume-averaged  simula¬ 
tions.  In  order  to  compare  different  simulation  techniques,  we 
studied  heat  transfer  behavior  such  as  the  volume  fraction  of 
melted  PCM,  average  temperature  at  different  locations,  tempera¬ 
ture  difference  between  micro-foam  and  PCM,  and  the  top-wall 
heat  fluxes. 

3.1.  Volume  fraction  of  melted  PCM 

Fig.  2  shows  the  profiles  of  melted  PCM  volume  fraction  with 
top  surface  temperature  fixed  at  334  K  and  368  K  respectively. 
Here  “1-T”  and  “2-T”  stand  for  one-  and  two-temperature  simula¬ 
tions.  “Arithmetic  Mean”  and  “Progelhoff’  refer  to  the  arithmetic 
mean  and  Progelhoff  models  for  the  effective  thermal  conductivity. 
First,  we  can  see  that  one-  and  two-temperature  simulations  with 
the  arithmetic  mean  predict  higher  melting  rate  than  those  from 
Progelhoff  model  and  the  DNS.  This  is  because  the  arithmetic  mean 
generates  higher  effective  thermal  conductivity,  thus  enables 
system  absorbing  more  heat  to  melt  PCM  quickly.  Second, 


Fig.  2.  Profiles  of  melted  PCM  volume  fraction:  (a )Ttop  =  334  K;  (b)  Ttop  =  368  I<. 
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one-temperature  and  two-temperature  simulations  accurately 
predict  volume  fraction  of  melted  PCM,  when  the  Progelhoff  model 
is  used  for  the  effective  thermal  conductivity  keff.  Hence,  modeling 
of  keff  play’s  a  key  role  in  the  volume-averaged  simulation.  As  pre¬ 
scribed  by  Progelhoff  model,  k^nonlinearly  depends  on  the  poros¬ 
ity,  and  is  dominated  by  the  thermal  conductivity  of  micro-foam. 
As  the  enforced  top  surface  temperature  is  raised  (Fig.  2b),  the 
temperature  gradient  at  the  top  surface  increases  and  induces  an 
earlier  melting  of  PCM. 

3.2.  Average  temperature  at  different  Y-cross  sections 

Fig.  3  shows  the  time  history  of  area-weighted  average  temper¬ 
ature  at  two  cross  sections  in  the  foam,  with  the  top-surface  tem¬ 
perature  fixed  at  334  K.  At  a  cross  section  close  to  the  top  surface 
(Fig.  3a),  the  temperature  simulated  from  volume-averaged  meth¬ 
od  and  Progelhoff  model  reasonably  fit  the  results  from  DNS.  Arith¬ 
metic  model  of  effective  thermal  conductivity,  however,  seemingly 
over-estimates  the  temperature  change  by  12%  in  the  phase 
change  region  compared  to  DNS. 

At  the  bottom  of  the  system  (Fig.  3b)  this  over-estimation  gets 
worse.  When  the  effective  thermal  conductivity  is  calculated  from 
arithmetic  mean,  the  volume-averaged  simulation  predicts  earlier 
starting  time  of  phase  change  compared  to  direct  simulation 


Fig.  3.  Average  temperature  profiles  at  (a)  cross  section  Y=  1.92  cm,  and  (b)  bottom 
plane  Y  =  0.0  cm.  The  top-surface  temperature  Ttop  =  334  K. 


(Fig.  3b),  indicating  the  over-estimation  of  the  effective  thermal 
conductivity  of  PCM/foam  system  by  arithmetic  mean.  When  the 
Progelhoff  model  is  applied  to  calculate  keff,  the  temperature  vari¬ 
ation  is  well  captured,  no  matter  whether  one-temperature  or  two- 
temperature  simulation  is  conducted.  It  is  noteworthy  that  there  is 
a  temperature  plateau  at  the  structure  bottom  (Fig.  3b),  in  which 
heat  transfer  is  primarily  dominated  by  phase  change.  In  contrast, 
the  average  temperature  at  the  cross  section  near  the  top  surface 
continuously  increases,  even  during  the  phase  change  process 
(Fig.  3a),  indicating  the  fast  rising  temperature  of  foam  structure 
at  this  cross  section  due  to  its  high  conductivity.  As  the  cross  sec¬ 
tion  moves  toward  the  bottom,  the  heat  conduction  through  foam 
structure  slows  down,  and  the  phase  change  process  is  more 
pronounced. 

3.3.  Temperature  difference  between  micro-foam  and  PCM 

In  two-  temperature  and  direct  numerical  simulations,  the 
effective  thermal  conductivity  and  the  interstitial  heat  transfer 
are  two  competing  factors  affecting  the  thermal  process  between 
micro-foam  and  PCM.  Larger  differences  in  values  of  effective  ther¬ 
mal  conductivities  between  micro-foam  and  PCM  induces  larger 
temperature  difference  between  them;  while  a  stronger  interstitial 
heat  transfer  makes  the  temperature  jump  smaller.  Comparison  of 
temperature  difference  cross  the  micro-foam  and  PCM  interface 
ST  =  TS-Tf  will  help  us  better  evaluate  the  two-temperature 
simulation. 

The  temperature  differences  at  two  cross  sections  (Y  =  1.92  cm 
and  0.0  cm)  from  two-temperature  and  direct  simulations  are 
shown  in  Fig.  4.  The  onset  of  phase  change  (melting)  is  indicated 
by  the  peak,  which  separates  pre-  and  post-melting  stages.  With 
results  from  direct  simulations  as  the  reference,  we  can  see  that 
the  two-temperature  simulation  using  arithmetic  mean  model  pre¬ 
dicts  earlier  melting  stages  in  both  Y-locations,  and  the  peak  mag¬ 
nitude  is  dramatically  over-estimated  at  the  cross  section  close  to 
the  top  surface  (Y=  1.92  cm).  However,  two-temperature  simula¬ 
tions  using  Progelhoff  model  predicts  acceptable  time  onset  of 
melting  process  and  the  temperature  difference.  With  the  cross 
section  further  away  from  the  top-wall  (e.g.,  at  Y  =  0.0  cm),  the  dis¬ 
crepancy  in  the  melting  onset  is  more  obvious. 

Temperature  difference  under  a  higher  top-wall  temperature 
( Ttop  =  368  K)  is  presented  in  Fig.  5.  The  melting  onset  is  advanced 
as  the  top-wall  temperature  is  increased.  We  can  see  that 


Fig.  4.  Comparison  of  temperature  difference  between  micro-foam  and  PCM 
(Top  =  334  K). 
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Fig.  5.  Comparison  of  temperature  difference  between  micro-foam  and  PCM 
(Ttop  =  368  I<). 


two-temperature  simulation  with  Progelhoff  model  predicts  nearly 
the  same  starting  time  of  phase  change  as  that  from  direct  simula¬ 
tion  in  both  cross  sections. 


3.4.  Heat  flux  at  top  surface 

We  further  calculated  the  area-weighted  average  heat  flux  at  top 
surface  using  various  simulation  methodologies.  In  one-tempera¬ 
ture  simulation,  only  the  overall  top-wall  heat  flux  can  be  calculated. 
In  two-temperature  simulation  and  DNS,  we  can  calculate  the  top- 
wall  heat  flux  for  each  component,  i.e.,  the  amount  of  heat  flows 
through  the  PCM  and  through  micro-foam  respectively.  The  overall 
heat  flux  at  the  top  surface  is  just  the  sum  of  heat  flux  from  PCM  and 
the  foam,  i.e.,  q"op  =  sq'ftop  +  (1  -  s)q's'[op,  where  s  is  the  area  ratio  of 
PCM  surface  area  to  the  total  on  the  top  surface.  The  heat  flux  of 
the  foam  is  defined  as  q''top  =  WVe^top-^p),  -m  rs  cd;@top  is  the 
foam  temperature  at  the  centroid  of  a  cell  near  the  top  surface,  and 
rd  is  the  distance  from  the  cell  centroid  to  the  top  surface.  In  the  cur¬ 
rent  study,  both  the  overall  and  individual  heat  fluxes  at  the  top  sur¬ 
face  are  normalized  by  the  upper  limit  heat  flux,  which  is  defined  as 
Qmax  =  H{ypfL  +  {ypfCps  +  (1  -  y)pscp,s)(Ttop  -  Tint)}.  H  is  the  dis- 
tance  from  the  top  to  the  bottom  surface,  and  is  the  cube  length  in 
our  case. 
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Fig.  6.  Heat  flux  at  top  surface  (a)  overall  (b)  individual  component  (Ttop  =  334  K). 
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Fig.  6(a)  shows  the  normalized  overall  heat  flux  at  top  surface. 
Volume-averaged  simulations  predict  the  heat  flux  reasonably 
well.  However,  the  effective  thermal  conductivity  calculated  from 
Progelhoff  model  generates  more  accurate  results  than  from  arith¬ 
metic  mean  model  (Fig.  6(a)  insert).  The  heat  flux  of  individual 
material  components  is  presented  in  Fig.  6(b).  From  the  insert,  it 
can  be  seen  that  the  heat  flux  from  the  micro-foam  is  much  larger 
than  that  from  PCM,  since  the  thermal  conductivity  of  the  micro¬ 
foam  is  three  orders  of  magnitude  higher  than  that  of  PCM. 
(Table  1). 

We  also  calculated  the  overall  and  individual  heat  fluxes  on  the 
top  surface,  at  which  the  temperature  Ttop  is  raised  to  368  K  (Fig.  7). 
The  flux  is  increased  at  the  top  surface,  contributed  principally  by 
the  micro-foam.  It  is  worth  mentioning  that,  unlike  the  tempera¬ 
ture  at  various  cross  sections,  the  heat  flux  gradually  and  slowly 
approaches  a  steady  value.  This  can  be  ascribed  to  the  varying 
effective  thermal  conductivity  and  temperature  gradient  in  the 
system. 

4.  Conclusions 

In  this  work,  we  numerically  studied  the  heat  transfer  behavior 
of  a  phase  change  material  in  an  open-cell  micro-foam  using  both 
volume-averaged  and  direct  simulations.  Heat  transfer  process  was 
simulated  for  a  spherical-pore  embedded  micro-foam,  and  its  top 
surface  was  kept  isothermal.  By  examining  the  volume  fraction 
of  melted  PCM,  average  temperature  at  different  locations,  temper¬ 
ature  jump  across  micro-foam  and  PCM  interface,  and  average  heat 
flux  at  the  top  surface,  we  find  that  volume-averaged  simulations 
can  reasonably  predict  the  heat  transfer  process  in  the  micro-foam 
with  complex  geometry,  as  long  as  accurate  thermo-physical 
parameters  such  as  effective  thermal  conductivity  are  provided. 
In  the  present  study,  the  effective  thermal  conductivity  of  the 
PCM/micro-foam  system  was  obtained  from  direct  numerical  sim¬ 
ulation  of  a  unit  cell.  The  simulations  results  supported  the  Progel¬ 
hoff  model  of  effective  thermal  conductivity  where  the  effective 
thermal  conductivity  can  be  expressed  as  a  power  law  of  porosity 
and  a  power  factor  of  1 .3  was  calculated  for  the  model  structure. 
We  found  that  the  heat  flux  near  the  enforced  isothermal  boundary 
is  dominated  by  the  micro-foam  due  to  its  high  thermal  conductiv¬ 
ity.  The  phase  change  heat  transfer  becomes  gradually  more  pro¬ 
nounced  in  locations  away  from  the  enforced  high  temperature 
boundary. 
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